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Abstract. Let A be a totally unimodular list of vectors in some lattice. Let 
Bx be the box spline defined by X. Its support is the zonotope Z(X). We 
show that any real-valued function defined on the set of lattice points in the 
interior of Z(X) can be extended to a function on Z(X) of the form p(D)Bx 
in a unique way, where p(D) is a differential operator that is contained in the 
so-called internal P-space. This was conjectured by Olga Holtz and Amos 
Ron. We also point out connections between this interpolation problem and 
matroid theory, including a deletion-contraction decomposition. 



1. Introduction 

Given a set = {u\, . . . , u^} of k distinct points on the real line and a function 
/ : — > M, it is well-known that there exists a unique polynomial pf in the space of 
univariate polynomials of degree at most k — 1 s. t. pf(ui) = f(v,i) for i = 1, . . . , k. 

If is contained in M. d for an integer d > 2. the situation becomes more difficult. 
Not all of the properties of the univariate case can be preserved simultaneously. 
The minimal number me s. t. for every / : — > M there exists a polynomial 
Pf E R[xi, . . . , Xd) of total degree at most me that satisfies Pf{ui) = /(u;) depends 
on the geometric configuration of the points in 0. Furthermore, the interpolating 
polynomial pj of degree at most me is in general not uniquely determined. This 
is only possible if the dimension of the space of polynomials of degree at most me 
happens to be equal to k. 

Uniqueness is possible if we choose the interpolating polynomials from a special 
space. Carl de Boor and Amos Ron introduced the least solution to the polynomial 
interpolation problem. For an arbitrary finite point set C M. d , they construct a 
space of multivariate polynomials 11(0) that has dimension |0| and that contains 
a unique polynomial interpolating polynomial pf for every function / : — > K 

H21CL3]. 

In this paper, we construct a space that contains unique interpolating functions 
for the special case where is the set of lattice points in the interior of a zonotope. 
The space is of a very special nature: it is obtained by applying certain differential 
operators to the box spline. This is interesting because it connects various algebraic 
and combinatorial structures with interpolation and approximation theory. 
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More information on multivariate polynomial interpolation can be found in the 
survey paper |17j . 

In this paper, we use the following setup: U denotes a <i-dimensioiial real vector 
space and A C U a lattice. Let X = (x\, . . . , xn) Q A be a finite list of vectors 
that spans U. We assume that X is totally unimodular with respect to A, i. e. 
every basis for U that can be selected from X is also a lattice basis. The symmetric 
algebra over U is denoted by Sym([/). We fix a basis Si, . . . , Sd for the lattice. This 
makes it possible to identify A with Z d , U with M d , Sym(f ) with the polynomial 
ring K[si, . . . , Sd], and X with a (d x N) matrix. Then X is totally unimodular if 
and only if every non-singular square submatrix of this matrix has determinant 1 
or —1. A base- free setup is however more convenient when working with quotient 
vector spaces. 

The zonotope Z{X) is defined as 

Z{X).= \^X i x i :Q<X i <l^. (1) 

We denote its set of interior lattice points by Z_(X) := int(Z(X)) n A. The box 
spline Bx ■ U — > K. is a piecewise polynomial function that is supported on the 
zonotope Z(X). It is defined by 

B x {u) := 1 vol w _ d { (Ai, . . . , \s) G [0, 1] N : ]T \x % = u) . (2) 

For examples, see Figure [1] and Example [TDJ A good reference for box splines and 
their applications in approximation theory is . Our terminology is closer to |14l 
Chapter 7], where splines are studied from an algebraic point of view. 

A vector u € U defines a linear form p x e Sym(£7). For a sublist Y C X, we 
define py := Y\ yeY Py For example, if Y = ((1,0), (1,2)), then py = s\ + 2s\S2- 
Now we define the 

central V-space V(X) := span{p y : rk(AT \ Y) = rk(X)} (3) 
and the internal V-space V-(X) := f] V{X \ x). (4) 

The space V- (X) was introduced in |18] where it was also shown that the dimension 
of this space is equal to |iJ_(AT)|. The space V(X) first appeared in approximation 
theory [T| fTUj. [TB]. Later, spaces of this type and generalisations were also studied 
by authors in other fields, e.g. [H H QH HOI HU [24] . 

We will let the elements of V-{X) act as differential operators on the box spline. 
For p £ V-(X) C Sym(J7) = M[si, . . . , s r ], we write p(D) to denote the differential 
operator obtained from p by replacing the variable Sj by 

The following proposition ensures that the box spline is sufficiently smooth so 
that the derivatives that appear in the Main Theorem actually exist. 

Proposition 1. Let XCACU = M. d bea list of vectors that is totally unimodular 
and let p G V-(X). Then p(D)Bx is a continuous function. 

Now we are ready to state the Main Theorem. It was conjectured by Olga Holtz 
and Amos Ron |18l Conjecture 1.8]. 
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FIGURE 1. A very simple two-dimensional example. Here, X = 
((1,0), (0, 1), (1, 1)), T-(X) = K, and \Z-(X)\ = 1. 

Theorem 2 (Main Theorem). Let XCACU = R d bea list of vectors that is 
totally unimodular. Let f be a real valued function on Z-{X), the set of interior 
lattice points of the zonotope defined by X . 

Then there exists a unique polynomial p £ V—(X) C R[si, . . . , s^], s.t. p(D)B\ 
equals f on Z-{X). 

Here, p(D) denotes the differential operator obtained from p by replacing the 
variable Sj by t^- and Bx denotes to the box spline defined by X . 

Remark 3. Total unimodularity of the list X is a crucial requirement in Theorem [5] 
Namely, the dimension of V-(X) and |-2L(X)| agree if and only if X is totally 
unimodular. Note that if one vector in X is multiplied by an integer A > 2, 
|,2L(X)| increases while V-(X) stays the same. 

Total unimodularity also enables us to make a simple deletion-contraction proof: 
it implies that A/a; is a lattice for all x G X. In general, quotients of lattices may 
contain torsion elements. 

Remark 4. We have mentioned above that dim("P-(A)) = |Z_(X)| holds. This 
is a consequence of a deep connection between the spaces V_(X) and V(X) and 
matroid theory. The Hilbert series of these two spaces are evaluations of the Tutte 
polynomial of the matroid defined by X [5] . One can deduce that the Hilbert series 
of the internal space is equal to the ^.-polynomial of the broken-circuit complex [5] of 
the matroid M*(X) that is dual to the matroid defined by X and the Hilbert series 
of the central space equals the /i-polynomial of the matroid complex of M*(X). The 
Ehrhart polynomial of a zonotope that is defined by a totally unimodular matrix is 
also an evaluation of the Tutte polynomial (see e. g. |25|). In summary, for a totally 
unimodular matrix X 

dim-p_(X) = \Z-(X)\ = Zx(0, 1) and <timV(X) = vol(Z(A)) = 1 x (l, 1) (5) 

holds, where 1x denotes the Tutte polynomial of the matroid defined by X. 

It is also interesting to know that the Ehrhart polynomial of an arbitrary zono- 
tope defined by an integer matrix is an evaluation of the arithmetic Tutte polyno- 
mial [61E]. 

Organisation of the article. In Section[5]we will discuss some basic properties of 
splines. We will prove the Main Theorem in the one-dimensional case in Section [3] 
In Section @] we will recall the wall-crossing formula for splines and employ it to 
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prove Proposition [TJ In Section [5] we will define deletion and contraction and prove 
two lemmas that will be used in Section [6] in the proof of the Main Theorem. 

2. Splines 

In this section we will introduce the multivariate spline and discuss some basic 
properties of splines. Proofs of the results that we mention here can be found in 
[Ffl Chapter 7] and some also in 

If the convex hull of the vectors in X does not contain 0, we define the multivari- 
ate spline (or truncated power) T x ■ U — > K by 

1 N 
T x (u) := = voljy_d{(Ai, . . . , Ajy) G R> : XjXj — u}. (6) 



^dct{XX T ) 



i=l 



The support of T x is the cone conc(X) := ^i x i '■ \ ^ o|. 

Sometimes it is useful to think of the two splines Bx and T x as distributions. 
In particular, one can then define the splines for lists X C U that do not span U. 

Remark 5. Let X C {/ ~ l r be a finite list of vectors. The multivariate spline Tx 
and the box spline Bx are distributions that are characterised by the formulae 

ip(u)B x (u)du= / ••• / ip iy~] XiXi dAi • • • d\ N (7) 



poo poo / \ 

/ " / ^ X] ' K%xi dA i ■ ■ ■ d ^ 

Jo Jo V i=1 j 



and J ip(u)Tx(u) du = j tp | } XjXj ) dAi ■ ■ ■ dAjy. 

where tp denotes a test function. 

Remark 6. Convolutions of splines are again splines. In particular, 

T x = T X1 * • • • * T XN and B x = B Xl * ■ ■ ■ * B Xn . (9) 

For x G X, differentiation of the two splines in direction x is particularly easy: 

D X T X = T X \ X (10) 
and D X B X = V X B X \ X ■= B x \ x - B X \ X (- - x). (11) 

Remark 7. For a basis C C U, 

D _ XZ(C) , T Xcone(C) , ^ 

Bc ~ c ~ Wicyy (12) 

where \A ■ U — > {0, 1} denotes the indicator function of the set A C U. In 
conjunction with ([9]), (|12[) provides a simple recursive method to calculate the 
splines. 

Remark 8. The box spline can easily be obtained from the multivariate spline. 
Namely, 

Bx{u) = (« " «s) , (13) 

sex 

where a s := J2a.es a - 
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FIGURE 2. The cardinal i?-splines B 2 , B 3 , and B 4 



3. Cardinal B-splines 

In this section we will prove Theorem [2] in the one-dimensional case. This will 
be the base case for the inductive proof of the Main Theorem in Section [6] 

Let Xn := (1, . . . , 1) C Z C K 1 . WLOG every totally unimodular list of vectors 

N times 

in IR 1 can be written in this way. 

One can easily calculate the corresponding box splines (cf. Remark [7]): 



B x N+1 (u) 



N+l 



3=0 



N\ 



where (u — j) 1 ? := max(u — j,0) N . The functions Bx N+1 are called cardinal B- 
splines in the literature (e. g. [S]). 

Note that Z_{X N+1 ) = {1,2,..., N}, 



Vx N+1 = span{l, 



r,N 



..,s 



}, and V-(Xn+i) = span{l, s, 



,,s 



,N-1 



}■ 



Hence, in the one-dimensional case, Theorem [5] is equivalent to the following pro- 
position. 

Proposition 9. Let N € N. For every function f : {1, . . . , N} — > M, there exist 
uniquely determined numbers \±, . . . , Xpf 6l s.t. 

N 

J2 KDl^Bx^ (j) = f(j) for j = 1, ... , N. (15) 

i=i 

Before proving this proposition, we give a few simple examples (see also Figure[2]). 
Example 10. 

B X2 (s) = s-2(s-l) + + (s-2) + 



Bx 3 (s) 



3(s 



2)1 



B Xi (s) =q( s3 - 4 ( s - !)+ + 6 ( s - 2 )+ 
The matrices M N are defined in (TlUl) below. 



(s - 3) 2 + ) 

4( S -3)^ + ( S -4) 3 + ) 



(16) 
(17) 

(18) 



M z = 1 



M 3 = 



_L 1 

2 2 



M 4 



'14 1 
6 6 6 
1 n _X 
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Proof of Proposition^ For N G N, we consider the matrix (TV x iV)-inatrix M N 
whose entries are given by 

m$ = D i - l B Xlf+1 (j). (19) 

The proposition is equivalent to AI N having full rank. The matrix M 2 = (1) 
obviously has full rank. Let us proceed by induction. By (fTTj) . D x Bx N+1 = 
V x D t ~ 1 Bx N - Thus, the matrices satisfy the following recursion: 

"»y = mf-i]j - m ?-i] 3 -i for * = !. • • • . N ~ 1. 3 = 1, • • • , N , and N > 2. (20) 

To simplify notation, we set m^ _1 = = 0. Let Vk, ■ • ■ , ujv denote the columns 

of M N . By induction, they are linearly independent. The columns of M N+1 are 
(ai, Vi — vq), • ■ • , (o-n + 1, vn+i — vn) with ay := Bx N+1 (j). We will now show that 
these vectors are linearly independent as well. Let Ai . . . , Ajv+i G R s. t. 

JV+l 

E x i°i = ( 21 ) 
i=i 

and X\(vi — vq) + ^2(^2 — «i) + . . . + Xn(vn — vn~i) — Xn+iVn = 0. The latter 
equation implies that all Ai are equal. We conclude that they must all be zero 
because of ([2"TT) and the fact that the oij are positive. □ 

4. Smoothness and wall-crossing 

The goal of this section is to prove Proposition [T] Before doing this, we mention 
some results on the structure of the multivariate spline Tx that are used in the 
proof. The Wall-Crossing Theorem describes the behaviour of Tx when we pass 
from one region of polynomiality to another. 

Definition 11. A tope is a connected component of the complement of 

Ux := {span(F) :Y CX, rk(Y) = rk(AT) -1}CU (22) 

The following theorem is a consequence of Lemma 3.3 and Proposition 3.7 in 
[Hi- 
Theorem 12. Let X C U = R d be a list of vectors N that spans U and whose 
convex hull does not contain 0. 

Then Tx agrees with a homogeneous polynomial f T of degree N — d on every 
tope t. 

Given a hyperplane H and a tope t which does not intersect H (but its closure 
may do so), we partition X \ H into two sets A T H and B T H . The set A T H contains 
the vectors that lie on the same side of H as r and B T H contains the vectors that 
lie on the other side. Note that the convex hull of {A T H , —BJj) does not contain 0. 
Hence, we can define the multivariate spline 

T x \h ~ (-l) lBJll T {A r H ,- BW . (23) 

Now we are ready to state the wall-crossing formula as in [15l Theorem 4.10]. 
Related results are in 151 1221. 



Theorem 13 (Wall-crossing for multivariate splines). Let n and T2 be two topes 
whose closures have an r — 1 dimensional intersection t\2 that spans a hyperplane 
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H. Then there exists a uniquely determined distribution f Tl2 that is supported on 
H s. t. the difference of the local pieces of Tx in t\ and ti is equal to the polynomial 

T£ - T? = (T X \ H T-^ H ) * r* . (24) 

Proof of Proposition^ We will show that p{D)Tx is continuous. By (fiU)) . this 
implies that p(D)Bx is continuous. We may always assume that is not contained 
in the convex hull of X: deleting zeroes from X changes neither P>x nor Z(X). In 
addition, one can always multiply a few vectors in X by —1 s.t. all vectors lie on 
one side of some hyperplane. This is equivalent to a translation of both, Z(X) and 
B x - 

Let u G U. If M G U \ Hx, there is nothing to prove: by Theorem [T2l Tx is 
polynomial in a neighbourhood of u and hence smooth. If u G "Hxi u is contained 
in the closure of at least two topes. We have to show that the derivatives of the 
polynomial pieces in the topes agree on u. This can be done using the wall-crossing 
formula. 

It is sufficient to prove that for two topes t\ and T2 that have an (r — 1)- 
dimensional intersection na, p(D)(T x 1 — T x 2 ) vanishes on T\2- 

Fix a vector x G X \ H. By definition, V-(X) C V(X \ x). This implies that 
p can be written as a linear combination of polynomials py where Y C X \ x and 
X \(Y U x) has rank d. Hence, X \ (H U Y) contains at least two vectors. 

By Theorem [H 

Dy(T x \ h - T-$) * r* = (T x \ { huy) ~ T~{l HUY) ) * DynH f Tl2 ■ (25) 
This polynomial is the convolution of a distribution supported on H with the dis- 
tribution (Tjcwjjuy") — ^x\(Huy))' ^ mce X\(H\JY) contains at least two elements, 
this polynomial vanishes on H. This finishes our proof. □ 

Remark 14. Holtz and Ron conjectured that V-(X) is spanned by polynomials py 
where Y runs over all sublists ofXs.t.X\(FU x) has full rank for all x G X |18l 
Conjecture 6.1]. By formula (fTl]) . this would have implied Proposition [T] However, 
this conjecture has recently been disproved [3]. 

5. Deletion and contraction 

In this section we will introduce the operations deletion and contraction which 
will be used in the proof of Theorem [5] in the next section. We will also prove two 
lemmas about deletion and contraction for box splines and zonotopes. 

Let x G X . We call the list X \ x the deletion of x. The image of X \ x under 
the canonical projection ir x : U — > [//span(x) =: U/x is called the contraction of 
x. It is denoted by X/x. 

The projection tt x induces a map Sym(7r :E ) : Sym([7) — > Sym(U/x). If we identify 
Sym({7) with the polynomial ring R[si, . . . , s r ] and x = s r , then Sym(7r a; ) is the 
map from K[si, . . . , s r ] to R[si, . . . , s r _i] that sends s r to zero and s\, . . . , s r _i to 
themselves. The space V(Xjx) is contained in the symmetric algebra Sym.(U/x). 

Since X is totally unimodular, A/a; C (7/j: is a lattice for every i£l and X/x 
is totally unimodular with respect to this lattice. 

Lemma 15. Let x G X , u G U , and u = u + span(a;) the coset of u in X/x. Then 
B x/X (u)= f B x \ x (u + tx)At = Y,Bx(u + \x). (26) 



8 



MATTHIAS LENZ 



Proof. Let <p : [//i->lbea test function and let ijj : U — >• M be a test function 
s. t. tp(u) = L ^(u + rx) dr. Note that a distribution T onU that is constant on all 
cosets of span(x) can be identified with a distribution T on U/x via T(tp) = T(tp). 
This is how (|26[) can be understood as an equality of distributions. We may assume 
that x = xn- Then 



/ ip(u)B x / x (u) du = / •••/ ip(j>~] XjXj) dXi ■ ■ ■ dAjy-i 
Ju/x Jo Jo , 



(27) 



^(T^AjXj +rx)drdAi ••■ dAjv-i (28) 



i 

tp(u) I B x \ x (u + Tx)dTdu. (29) 
Jr 

This proves the first equality. For the second equality, note that 

/ B x \ x (u + tx) dr = / B x \x(u + Xx — tx) dr = B x {u + Ax). □ 

AGZ AGZ 

Remark 16. Lemma[Tn]is a statement on semi-discrete and continuous convolutions 
with the box spline. A related result is in [23] . 

The following lemma yields a deletion-contraction formula for the interior points 
of the zonotope. See Figure [3] for an illustration. 

Lemma 17. The following map is a bijection: 

Z-(X) \Z-(X \x) -> Z-(X/x) (30) 

z i y z . (31) 

Proof. It is obvious that z is contained in Z-(X/x). Using the fact that Z_(X)| 
is an evaluation of the Tutte polynomial (formula ([5])) and the deletion-contraction 
formula for the Tutte polynomial, one can easily establish that the domain and the 
range of the map have the same cardinality. 

Hence it is sufficient to show that the map is injective. Let us prove this. First 
note that z G Z-(X) is contained Z_(X \ x) if and only if z + x G Z-(X). Let 
zi,Z2 G Z-(X) \ Z^(X \ x) s. t. z\ = z<i- This implies that there is a A G R s.t. 
z\ = Z2 + Ax. Because of the total unimodularity, A must be an integer. WLOG A 
is non-negative. By convexity, Z\, . . . , z\ + x, . . . , Z\ + Xx are contained in Z-(X). 
By the observation at the beginning of this paragraph, z\ + x is not contained in 
Z-(X). This implies A = 0. Hence the map is injective. □ 

6. Exact sequences 

In this section we will prove the Main Theorem. We start with a simple obser- 
vation. 

Remark 18. If X contains a coloop, i. e. an element x s.t. rk(X \ x) < rk(X), then 
Z..{X) = and V-(X) = {0}. Hence, Theorem |2] is trivially satisfied. 

We will consider the set S(X) := {/ : A — > K : supp(/) C Z_(A)} and the map 

lx:V-(X)^S(X) 

p^[A9 z^p{D)B x (z)] 



INTERPOLATION, BOX SPLINES, AND LATTICE POINTS IN ZONOTOPES 




Figure 3. Deletion and contraction for a zonotope and a function 
defined on the interior lattice points of the zonotope. 



Proposition 19. Let d > 2 and let A C U = R d be a lattice. Let X C A be a finite 
list of vectors that spans U and that is totally unimodular with respect to A. Let 
x G X be a non-zero element s. t. rk(AT \ x) = rk(X). 

Then the following diagram of real vector spaces is commutative, the rows are 
exact and the vertical maps are isomorphisms: 







■ r^(x\x)-^r-(x) s -^K(x/x) — -o (33) 



7x\n 



V, 



7x/i 



^ E(X \ x) — =— >- S(X) E(X/x) ^ 

where V x (f )(z) := f(z)- f(z-x), (34) 
^x{f){z):= ^2 f(x) = ^2f(\x + z) for some z e z. (35) 

irgznA AeZ 

Proof. Commutativity of the left square: Let z G Z-{X) and let p £ V-{X \ x). 
By (HU, {p- Px )(D)B x {z) = V x { P (D)B x \ x )(z). Hence lx ° (-p x ) - V x o lx \ x . 

Commutativity of the right square: Let z S Z_(X/x), f G V-(X) and let z G U 
be a representative of z. Then 

Y J P(D)B x (\x + z)=p(D) f B XVx (u + Tx)dT = Sym(7r x )(p)(D)B x/x (z) 

because of Lemma ll5l and the fact that applying a differential operator to a function 
that is constant on a subspace is the same as applying the projection of the differ- 
ential operator to the projection of the function. Hence *y x / x o Sym(7r T ) = T, x o r y x . 

Exactness: Exactness of the first row was stated in [2] and proven in |20] . The 
proof relies on the fact that V-(X) can be written as the kernel of a power ideal. 
Exactness of the second row is easy to check taking into account Lemma 1171 

Isomorphisms: By induction over the number of non-zero elements in X, r y x \ x 
and "f x / x are isomorphisms. Then y x is a l so an isomorphism by the five lemma. 

Two base cases have to be considered: by deleting elements from X, it may 
happen that X eventually contains only coloops and zeroes. This case is trivial 
(cf. Remark [HI). 

By contracting elements from X. it may happen that X has rank 1. We have 
shown in Section [3] that V~(X) and S(AT) are isomorphic in this case. □ 

Proof of the Main Theorem. The theorem is equivalent to "f X being an isomorph- 
ism, which is part of Proposition [121 D 
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